Hawking radiation for a Proca field in D dimensions II: 
charged field in a brane charged black hole 
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We generalise our first analysis of the wave equation for a massive vector boson in the background 
of a D-dimensional Schwarzschild black hole p]], by adding charge both to the field and the black hole, 
on the 3+1 dimensional Standard Model brane. A detailed numerical study is performed to obtain 
the transmission factor for the coupled (as well as decoupled) system of equations describing the 
Proca field modes, varying the angular momentum number, mass, charge and space-time dimensions. 
A qualitatively new feature arising from the introduction of charge is the appearance of superradiant 
modes, which we investigate. We then compute the Hawking fluxes, and analyse the effect of the 
charge. In particular we observe an inverted charge splitting effect for small energies and for two or 
more extra dimensions. For neutral fields, we compare the emission of massive particles with spins 
up to one and also compare the Proca bulk-to-brane ratio of energy emission, showing that, as for 
a scalar field, most of the energy is emitted on the brane. 

PACS numbers: 04.50.-h, 04.50.Kd, 04.20.Jb 



I. INTRODUCTION 

The possible existence of additional space-like dimen- 
sions in nature has been discussed for a century, at least 
since the works by Nordstrom, Kaluza and Klein. Dur- 
ing the last four decades, moreover, the naturalness of 
extra dimensions within Supergravity and String theory 
made it a main stream topic within high energy theoret- 
ical physics. At the end of the last century, this research 
led to models that, aiming at solving the hierarchy prob- 
lem, predicted that the extra dimensions could be very 
large (or even infinite), as compared to the traditional 
Planck scale. Within such scenario, the true fundamen- 
tal Planck scale could be as low as the TeV scale [2- 5] so 
that strong gravity effects could be visible in realistic man 
made particle accelerators. In particular, miniature black 
holes would be produced in particle collisions and Hawk- 
ing radiation would be the main observable signature 

. This motivation led to an intensive study of Hawking 
radiation from higher-dimensional black holes [9I-I2TI . 

Though so far there is no experimental evidence for 
the existence of such extra dimensions or strong gravity 
events in the currently probed energy frontier at short 
distances, the bounds are not yet as strong as one might 
wish |23 - |27| . Thus, it is still of interest to understand 
how such events should reveal themselves, as to impose 
the best possible bounds [H, [2^, especially for future 
high energy experiments, where a semi-classical treat- 
ment of the Hawking evaporation becomes a better ap- 
proximation. Such an approximation holds in the highly 
trans-planckian regime, otherwise quantum gravity cor- 
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rections become important |30l-l3 

One of the Hawking radiation channels that has not 
been properly addressed in the literature is that of mas- 
sive vector bosons, both uncharged and charged, to de- 
scribe the emission of Z and particles of the Stan- 
dard Model (SM). The basic difficulty is that the Proca 
equations do not decouple even in a spherically symmet- 
ric black hole background. To bridge this gap, we have 
recently constructed a numerical strategy to investigate 
the transmission factors, as well as the Hawking radia- 
tion flux, for a neutral Proca field in a 13 dimensional 
Schwarzschild s pac e-time In the present paper we 
close this study j20l [2ll [ by considering the Proca field on 
a brane and including charge for both the background 
and the spin 1 field, as to model the charged vector 
bosons of the SM. 

A qualitatively new feature that occurs when both 
background and Proca field charges are included is the 
existence of superradiant modes. For a bosonic field 
mode with frequency cj and azimuthal quantum number 
TO, it is expected to observe superradiant scattering in the 
background of a rotating black hole with horizon angular 
velocity fin when uj < rnQn- These modes are amplified 
by extracting rotational energy and angular momentum 
from the black hole. A similar amplification effect can be 
observed for a charged bosonic field mode, with frequency 
id and charge q, when impinging in a charged black hole 
with electric potential at the horizon ^h, when uj < q^H- 
Such mode extracts Coulomb energy, as well as charge, 
from the charged black hole. In a rotating background, 
the Proca field equation variables are not known to sepa- 
rate, which presents an extra difficulty added to the non 
decoupling of the modes, making it difficult to study ex- 
actly the superradiance phenomenon - see 33, 34] for a 
recent study in the slow rotation approximation. Thus, 
as a second motivation to our work, the charged black 
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hole background, due to the spherical symmetry, yields a 
setup where superradiance of a massive spin 1 field can be 
studied without any approximation, albeit numerically. 
Such analysis will be performed herein. 

The structure of this paper is organised as follows. In 
Section ini we introduce the background geometry which 
will be explored in this paper. In Section IlIII we present 
the equation of motion for a charged Proca field on the 
charged brane; we discuss the boundary conditions and 
its conversion to a first order system adequate for nu- 
merical integration in Section llVI The numerical results 
for the transmission factor and the associated Hawking 
fluxes are presented in Section |V] and we close with a 
discussion and final remarks. To keep the main part of 
this paper compact and clear, some technical relations 
are left to an Appendix. 

We shall use the particle physics signature convention 
(H ) in this paper. 



II. PROCA FIELD AND BACKGROUND 



A. Massive W 

There are two sub-cases which have to be dealt with 
separately. Modes with i ^ 0, where £ is the usual an- 
gular momentum quantum number, are described by two 
coupled radial fields {ip,x)^ corresponding to two prop- 
agating deg rees of freedom and a decoupled transverse 
field T [3(|. The radial equations for the coupled fields 
are 







ar 





+ [2iVujr - ir (cor - qQ) V']x = ; 

V'r2— + {ur - qQf - V + 1) + Af^r^) 
+ [2iqQV - ir [uor - qQ) V']'ip ^ Q , 



(3) 



whereas the decoupled equation for the transverse mode 
is 



We consider a Proca field with mass M and charge 
q, on a charged brane, which describes the SM W parti- 
cle. The Lagrangian for such field in the SM of particle 
physics (after electroweak symmetry breaking in unitary 
gauge) is 



iqWlW.F^^" , (1) 



where W,, 



and we 



have included the coupling of the W field to the elec- 
tromagnetic field strength tensor — d^A^, — d„A^ 
as allowed by gauge invariance in the SM. The electro- 
magnetic potential is denoted as A^^. We consider the 
following background brane black hole geometry 

1 



ds^ = V{r)dt^ 



V{r) 

with metric function 

V{r) = 1 



dr^-r^(dr -t-sin^ 61 V) , (2) 



r2 



where fi and Q are the parameters related with black hole 
bulk mass and brane charge. For numerical convenience, 
we choose units such that the outer horizon radius is 
rn = 1, i.e. fi = 1 + at the outer horizon. 



III. SYSTEM OF RADIAL EQUATIONS 

Due to the spherical symmetry of the background, we 
can obtain the equations of motion, using the gauge- 
invariant method of [35], following the same procedure 
as we have done in [l| . The form of the equations is very 
similar to that obtained therein; we therefore present here 
only a summary of the various radial equations (after 
Fourier transforming with respect to time and factoring 
out the angular harmonics). 



T = 0. 

(4) 



Here £ starts from one, and uj is the frequency of the 
modes. For the special case £ — 0, there is a single mode, 
as expected, and the two other modes are non-dynamical 
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B. Massless 

A special case, whose modes must match two of the 
modes in the previous subsection when M — 0, is the 
exactly massless theory. This is an interesting case for 
various reasons. Firstly it provides a way to check the re- 
sults of solving the coupled system in the previous section 
in a special limit. Secondly, it must provide accurate re- 
sults for two of the degrees of freedom of the Proca field 
at high energies. Finally, the radial equations for this 
special case are all decoupled, so the definition of trans- 
mission factors is straightforward. Again this works as a 
numeric check of the results for the transmission factors 
of the coupled fields when Af ^ 0. 

In this special case, the £ = Q mode is non-dynamical, 
i.e. its wave equations are satisfied identically. For ^ 7^ 0, 
there is a transverse vector mode which is obtained by 
setting M = in Eq. and another mode described 
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by X which satisfies 



dr V dr 



_2qQVh_d_ 

{ujr — qQ) dr 



[ujr - qQf - l{i + l)V 



X = 0, (5) 



whereas f/' = irVdrxl {'^^ ~ qQ) non-dynamical. We 
would like to anticipate that despite the apparent difi^er- 
ent form for these two modes, x and T, we shall find 
they give the same numerical results for the transmission 
factor, in all dimensions. 

Since decoupled radial equations have been extensively 
studied in the literature, and a similar coupled system 
has been studied in detail in [l| , in the remainder we will 
only review some of the details of the numerical set up 
to solve the coupled system. 



IV. BOUNDARY CONDITIONS AND FIRST 
ORDER SYSTEM 

To determine the transmission factors, we need to in- 
tegrate the radial equations from the horizon to the far 
away region with ingoing boundary conditions. The stan- 
dard procedure is to find a series expansion of the solution 
near the horizon, which can be used to initialise the so- 
lution (we do so at r = 1.001). Focusing on the coupled 
system {V', x}: if we define y — r — 1, Eqs. ([3]) become 



^dy^ 



^(0:^ + -8(0-^ + ^(0 



dy 
d 
' dy 



i> + E{T)x^Q, (6) 
X + E{T)ij^{), (7) 



where the polynomials A, B, C, E, A, B, C, E are defined 
in the Appendix. Making use of Frobenius' method to 
expand ^ and Xj we insert the following expansions into 
Eqs. © and dZl) 

OO OO 



3=0 



(8) 

where the sign of p was chosen to impose an ingoing 
boundary condition, and we obtain the recurrence re- 
lations (|Aip for the coefficients /ij and Vj found in the 
Appendix. A general solution close to the horizon can be 
parametrised by two free coefficients and vi . 

Similarly, to understand the asymptotic behaviour of 
the waves at infinity we now expand ip and x as 



X = e''r^T.^^ (9) 

j=o ^ 



which after insertion into Eqs. ^ yield 
/3 = ±ik , p = ±iip , 



(10) 



where (p = (5„,2(w^ + fc2)(l + Q'^)/{2k) - qQco/k, and 
k = — M^. Thus one can show that asymptotically 
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(12) 



where $ = /cr + (^ log r and is defined in the Appendix, 
Eq. (|A2p . Thus, as expected, each field is a combination 
of ingoing and outgoing waves at infinity. Asymptoti- 
cally, the solution is parametrised by four independent 
coefficients {0^,05*^}, two for each independent mode in 
the coupled system. In the same way as in [l|, one can 
define a first order system of ODEs containing four radial 
functions {x^i V'^} which coincide with such coefficients 
at infinity, allowing for an easy extraction of the wave 
amplitudes. Our target system which will be solved nu- 
merically in the remainder is then 



dr 



XT 



dT 

dr 



(13) 



with — ('0+, ■(/)_, x+, X-)- The definition of the ma- 
trices X and T, and how they relate with (|6]) and Q are 
found in the Appendix. 



V. HAWKING FLUXES 

We shall now calculate the transmission factor for the 
coupled system as well as the Hawking fluxes generated 
from all the modes. In [l| we have found a conserved 
flux from the energy momentum tensor which was simple 
enough to extract the transmission factors. However, in 
the present case, the energy-momentum tensor is more 
difficult to deal with so we use instead the conserved 
electric current which is naturally defined for this charged 
field. One can show that such a current is given by 



c.c. 



(14) 



The radial flux at r fixed is obtained by integrating the 
a ~ r component on the sphere. We note that only 
the flrst term (denoted Jf) contributes, since the second 
term becomes a total derivative on the sphere. 

The contribution for the flux for the coupled fields at 
infinity, is then found by simplifying the radial compo- 
nent of Eq. (jl4p using the equations of motion, and in- 
serting the far away expansion at infinity 
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(15) 
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where y|(s = ±;i = 0,1) are hnear combinations of 
the asymptotic coefficients af given in the Appendix, 
Eq. (jA7|) . Using the reflection matrix (R) defined in [l|, 
we obtain 

JiLZupie - (y-)^ (1 -RtR)y- ^ (y-)tTy- , (16) 

where we have defined a (hermitian) transmission matrix 
T, which can be diagonalised such that the decoupled 
asymptotic fields are found. 

Following the same procedure, one can also calculate 
the electric current at the horizon 



^ I — couple 



1 



Lj ~ qQ 



(h-)^h-. 



(17) 



where the coefficients are linear combinations of the 
two independent coefficients (i = 0, 1), given in the Ap- 
pendix, Eqs. (jASp . Eq. ([T7]) shows the important point 
that the current can be positive or negative. This is ex- 
pected, because for the boson field, the electric coupling 
can trigger superradiance. 

Furthermore, from the conservation law of the electric 
current, one can find an alternative expression for the 
transmission matrix. Using the scattering matrix (S ) 
defined in [1], we find 



1 



uj - qQ 



(s— St 



— ot — \-l 



(18) 



Once we have obtained the transmission factors, the 
number and energy fluxes are given by 



d{N,E} _ 1 ^ 



dtduj 



AC 



where C labels the mode and Th is the Hawking temper- 
ature which, in our units, is 



H 



(n - 1) + (n - 3)Q^ 
47r 



(20) 



VI. NUMERICAL RESULTS 



We now present a selection of numerical results for 
the transmission factor and the corresponding Hawking 
fluxes. In order to integrate the decoupled and coupled 
radial equations, we wrote independent codes in math- 
EMATICA7 and in CH — h, finding agreement between the 
two codes. Using them we have generated a set of figures 
that we now describe. 

In Fig. [1] the transmission factors for different masses, 
space-time dimensions and charges are illustrated to ex- 
hibit the superradiance phenomenon. As mentioned in 
the introduction, superradiant amplification of a bosonic 
field in a charged and/or rotating black hole occurs since 
there is Coulomb and/or rotational energy that can be 
extracted without decreasing the black hole area. The 
general condition of superradiance is a; < rnQn + Q^h', 



when this condition is verified the transmission factor be- 
comes negative and the scattered mode is amplified. In 
order to make the results clearer and more readable, here 
we just show the two coupled modes with £ ^ 1, as an 
example. Other cases are qualitatively similar. In the 
left panel of Fig. [TJ we show the superradiance depen- 
dence on the field charge for n = 2 (top row) and n = 4 
(bottom row) . It is clear that superradiant amplification 
is enhanced with growing field charge except in the small 
energy regime; in this regime one observes, at least for 
the cases plotted, that the amplification decreases with 
increasing field charge (this effect is more noticeable for 
higher dimensions - bottom row). The middle panel of 
Fig. [T] shows the dependence on the background charge 
- the trend is the same as when the field charge is var- 
ied, i.e. a generic superradiance enhancement with the 
charge except for small energies. In the right panel of 
Fig [TJ we show the superradiance dependence on the 
field mass; a suppression effect with increasing mass is 
observed, which is a generic behaviour for the transmis- 
sion factor (independently of the spin of the field or of 
being in the superradiant regime) (20l . [2l| . 

After we obtain the transmission factors, the Hawking 
fiuxes (fT9|) can be calculated directly. In the remainder 
we have chosen to show the fiux of the number of parti- 
cles. The flux of energy has similar features and can be 
simply calculated through multiplying each point by w. 

In Fig [2] we illustrate the number fluxes dependence 
on the field charge. On the left plots, we have kept the 
background charge parameter fixed (and positive), with 
Q — 0.4 and n = 4, varying the field charge q from — 1 
to 1. The plots show that there is a region at low energy 
where the negative field charge flux is larger, whereas 
in the remaining spectrum positive charge emission is 
favoured. It is also clear that if we integrate over the 
curves the emission of positive charge is always favoured. 
This low energy behaviour where negative charges are 
favoured, only occurs for n > 3 as can be seen on the right 
panel, where the difference between the positive charge 
and negative charge flux spectrum is presented for various 
n. This inverted charge splitting effect was also observed 
for scalars and fermions [2l|, and it results from the in- 
terplay between the thermal factor and the transmission 
factor appearing in the expression for the number fluxes. 
Whereas the thermal factor always favours same charge 
emission, the transmission factor favours opposite charge 
emission and these factors dominate different parts of the 
spectrum. We have considered the same parameters as 
in [2l!|, Q = \q\ = 0.3, wherein scalars and fermions have 
been studied, to allow for an easy comparison. 

In Fig. [3l we present the number flux dependence 
on the background charge for different field charge in 
the same row and different space-time dimensions in 
the same column. Consider first the q — case 
(middle column); it shows that the fluxes are sup- 
pressed/maintained/enhanced with the increase of back- 
ground charge for n = 2/n = 3/n = A. To understand 
this behaviour observe, from the definition of Hawking 
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FIG. 1. Negative transmission factors that show the superradiant amphfication of the coupled modes with 1=1 for different 
parameters. The two coupled modes are represented with the same colour (solid and dashed lines). The top and bottom rows 
differ in the space-time dimension. (Left panel) Variation with the field charge q; (Middle panel) variation with the background 
charge Q\ (Right panel) variation with the field mass. 




FIG. 2. (Left panel) Variation of the number fluxes with field charge for fixed n and black hole parameters in the small field 
mass limit. (Right panel) Variation of the difference between positive charge and negative charge flux with n. 



temperature, Eq. (j20p . that the Hawking temperature 
is decreased/maintained/increased with increasing back- 
ground charge for n = 2/n = 3/n = 4, in these units. 
Since we are using horizon radius units, as we vary the 
charge parameter Q, we are actually varying the mass of 



the black hole as well as the charge while keeping = ^■ 
Nevertheless, it is easy to see that, up to a stretching of 
the horizontal axis, if we fix the black hole mass and vary 
the dimensionful charge, these conclusions for the varia- 
tion of the height of the curves do not change since the 
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dN 
dujdt 



0.008 - 
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FIG. 3. Number fluxes dependence on the background charge for different space-time dimensions and field charge, with fixed 
field mass M — 0.4. We vary the field charge in the same row and vary the space-time dimension in the same column. 



number flux is dimensionless [38|. For higher tempera- 
ture, one expects a larger flux of particles, which is in- 
deed the behaviour shown in the second column of Fig. |3l 
Turning on the field charge we observe a more involved 
behaviour. For n = 3, in which the Hawking temperature 
does not vary with Q, we see in the first/third column 
and for sufficiently large energies a monotonic suppres- 
sion/enhancement of the Hawking flux when the black 
hole has the opposite/same charge as the field. This 
is in agreement with the discussion of the left panel of 
Figure [2] For n = 2,4, varying Q one also varies the 
Hawking temperature and more complex patterns are 



observed. Another trend is that the number fluxes in- 
crease as the space-time dimension increases, which may 
be understood from the existence of more modes that 
contribute to the transmission factor. 



A. Mass effect and bulk/brane emission 

In Figure |4l we perform a comparison of the effect of in- 
troducing a mass term for the various spins which are rel- 
evant for the Standard Model brane degrees of freedom. 
We have used the data of ^] for scalars and fermions. 
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n = 4, g = Q = 




FIG. 4. Variation of the number fiuxes for neutral particles with different spins (scalar, fermion and Proca) on the brane. For 
each type of particles we consider three different masses (M = 0, 0.5 and f), which can be identified by the starting point of 
the curve. Note we have included the two helicities for the fermion field and all the three modes for the Proca field (we have 
used a small mass M = 0.01 for the latter, instead of M = 0). 



Note that for fermions we have multiplied the data by 
a factor of two to take into account the two helicities 
of the Dirac field, since here we are also considering all 
the three modes for the Proca field. The dashed curves 
are for increasingly larger field mass (we have used the 
cases M = 0, 0.5 and 1, as can be seen from the thresh- 
old points where the curves start). For D = n + 2 = 4 
(left panel), we observe a striking similarity between the 
Proca (s = 1) flux with the scalar (s = 0) flux for AI — 
(except at the high energy tail) , which is due to the dom- 
inance of the ^ = mode. We can see this feature is 
always true at small energies for larger n (middle and 
right panels); as we increase n, however, higher modes of 
the Proca field enhance the fiux as compared to scalars. 
The extra modes of the Proca field contributing at higher 
energy, also explain the fact that the mass suppression 
is not as large as for Dirac fermions or scalars, as we 
see from the dashed curves corresponding to M = 0.5 
for example. Contrasting with the low energy behaviour, 
where scalars and Proca field are dominated by the s- 
wave, whereas Dirac fermions are suppressed since they 
do not allow an s-wave, is the high energy behaviour. 
There we observe tails which are in the ratio 1:2:3 fol- 
lowing the number of degrees of freedom for the scalar, 
fermion and Proca fields respectively. This is in agree- 
ment with the fact that all transmission factors tend to 
one at high energy. 

In the remainder, we combine our results for the Proca 
field on the brane with those in (Tj for a Proca field in 
the bulk, to analyse the relative bulk-to-brane emissivity 
in the neutral case. 

The total energy rate (or number rate) emitted in 
Hawking radiation for a given field, is given by integrat- 
ing the fluxes of Eq. p9)) (or their counterparts in the 
bulk) over w. The bulk-to-brane energy emissivity ratio, 
for massive neutral Proca fields of different mass is shown 
in TableUas a function of n. We have used the data in [ij 





n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


M = 0.1 


1 


0.46 


0.38 


0.41 


0.53 


0.83 


M = 0.3 


1 


0.49 


0.40 


0.42 


0.55 


0.86 


M = 0.6 


1 


0.59 


0.47 


0.49 


0.62 


0.95 



TABLE I. Bulk-to-Brane relative energy emission rates for 
massive neutral vector fields for different mass M in terms of 
the space-time dimension D — 2 + n. 



to obtain the total emissivity for bulk fields, as well as the 
data presented here, for the brane emissivity. The entries 
of Table HI show clearly that the emission of energy into 
brane- localised Proca particles is dominant, for n larger 
than two, which is consistent with the argument that 
black holes radiate mainly on the brane For fixed 

mass, as n increases, the bulk-to-brane energy emission 
ratio initially decreases, reaching a minimum value for an 
intermediate n. However, if one increases n further, the 
bulk-to-brane energy emission ratio increases again. Fur- 
thermore, we have also observed that the bulk-to-brane 
energy emission ratio increases with the field mass. A 
similar behaviour was observed for scalar fields [l^, 113 ■ 
Finally, we found that the bulk-to-brane energy emission 
rate for the Proca field is larger than that for the scalar 
field for all n, being more noticeable for large n, say n = 6 
and n = 7 (note the different notation for D = 2 + n in 
our paper as compared with [lol l40j). 



VII. DISCUSSION AND FINAL REMARKS 

In this paper we have completed the analysis started 
in [H , by computing the transmission factor for a charged 
Proca field propagating in the background of a charged 
black hole on a brane. Furthermore, this study completes 
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the study of the effect of mass and charge for particles 
evaporating on the brane, for all spins relevant for the 
Standard Model in brane world scenarios [IS, HH . Since 
the Proca field equations for the various modes do not 
totally decouple, we have followed a numerical strategy, 
to determine the transmission factors as in 

One of the main, and novel features arising from con- 
sidering a charged Proca field on a charged background is 
the existence of negative transmission factors, a signature 
of superradiant scattering, which we have presented and 
described. Our results are consistent with the condition 
of superradiance as allowed by the area theorem, and we 
found that generically increasing the field or background 
charge amplifies the effect (though some inversions are 
possible at small energies). It is worth commenting that 
although superradiant scattering is observed for both 
rotating and charged black holes, rotating black holes 
exhibit superradiant instabilities against massive scalar 
fields, and in particular against the Proca field fssl. Is^j. 
whereas charged (non-rotating) black holes do not seem 
to exhibit analogous instabilities against charged, mas- 
sive bosonic fields. This can be checked for the Proca 
field using our numerical framework, by computing the 
frequencies of quasi-bound states, a study that can also 
teach us how long lived Proca hair around charged black 
holes can be. Thus, even though our setup offers the 
technical advantage of spherical symmetry for studying 
superradiant scattering of the Proca field, it cannot be 
used for studying superradiant instabilities, and in par- 
ticular their non-linear development. Given the lack of 
spherical symmetry in the rotating case, such study will 
require sophisticated numerical relativity methods of the 
sort used in [41-43]. 

An effect observed herein, is that similarly to scalars 
and fermions, there is an inverted charge splitting ef- 
fect at small energies for more than one extra space- 
time dimension, where particles with charge opposite to 
the black hole are emitted dominantly. Nevertheless, 
for larger energies, the normal splitting order, favour- 
ing the emission of same charge particles as to discharge 
the black hole is restored and overall (by integrating out 
the flux), this channel tends to discharge the black hole. 
This effect has been suggested to be another signature 
that could be found in the energy spectrum of charged 
fermions [2l|. Since fermions in the final state are also 



produced indirectly from the decay of vector bosons such 
as or the Z particles one may ask whether the effect 
survives. Here we have verified that the effect is present 
for W^, so at least for the case when the final state de- 
cay products are a charged lepton £^ and a neutrino, this 
contributions will certainly enhance the effect. 

In the neutral case, we have performed two compar- 
isons. First we have compared the effect of the mass and 
spin on the Hawking fluxes, for all spins in the Stan- 
dard model. Our main findings are that the Proca field 
spectrum departs from being similar to a scalar field in 
four dimensions and becomes increasingly dominant for 
larger number of dimensions, peaking and extending to- 



wards larger energies. This also means that the mass 
suppression effect is smaller for the Proca field than it is 
for scalars and fermions. Secondly, we have compared the 
bulk-to-brane emission ratio for Proca fields, confirming 
brane dominance in general, and the suppression of brane 
dominance with increasing mass. 

Our methods and results can be used to improve the 
modelling of black hole evaporation in TeV gravity sce- 
narios, in the black hole event generators psi . |29| that 
are in use at the ATLAS and CMS experiments to put 
bounds on extra dimensions in this channel [22h241 [2711 . 



Acknowledgements. M.W. and M.S. are funded by 
FCT through the grants SFRH/BD/51648/2011 and 
SFRH/BPD/69971/2010. The work in this paper is also 
supported by the grants PTDC/FIS/116625/2010 and 
NRHEP-295189-FP7-PEOPLE-2011-IRSES. 



Appendix A: Functions and matrices 

The functions that are used in the text are (where Kq 

iii+l)): 

2n+l 

A{r) = = r [r" - (1 + + QV"-^] ' , 

m=0 



B{r) ^ = 2 [r" - (1 



m=0 



2n+l 



C{r) = Y c^ny"^ ^ (^'- - iQfr^"'^ 



ni—0 



- {kI + MV2)r"-i [r" - (1 + Q2)r + gV"-^] , 
E{r) = Y ^^V"" = i^Q^""' [(1 + - 1)^ 

m=0 



2n 



A{r) = Y amy" = [r" - (1 + Q'y + QV"-^]' 



m— 



B{r) = , 



C{r) = Y Cmy™ = i^r - qQfr^^- 

m=0 



9 



2n-2 



E{t) = e^y^ = i(?Qr"-i [{n - 3)(1 + Q2)+ 
+2r"-i] - iwr"-i [(1 + Q2)(„ _ i)^ _ 2QV""2] 

The recurrence relations are 



Ml = 



[azpjp - 1) + 62P + ci + ei] i/Q + epi/i 
a2/5(p + 1) + Co 



with 



a2(p + j)(p + j - 1) + co X 

X [a2(p + j)(p + j - 1) + Co] - eoeo , 

j 

Y [(am+2 (p + i - to) (p + i - m - 1) 

m— 1 

m+1 

(p + j - m) + Cm) Mi 



_ a2(M + i)(p + i - 1) + go ^ 



a2(p + j)(M + j - 1) + Co ? Co 



[(a™+2(M + i - to)(p + j - TO - 1) + c™) 



m— 1 



/.-^/., (Al) 



J 



The coefhcients used in the asymptotic expansion in the text are 



2u 



q^Q^]VP qQuJ 



2fc 



2k 



(A2) 



The relation between the new first order radial functions used in (|13p . and the 4- vector V-^ = (V', drV': X^ drx) for 
the original fields and derivatives is found from Eqs. pT|) . (IT^ and their derivatives. The corresponding r-dependent 
matrix transformation T is defined 



and its form is 



£1 




r 


r 






r 


r 















V = T* 



(A3) 



[ik + f ] e^* 



-1$ 



(A4) 



On another hand, the original system ([3]), can be written in a first order form 

dV 



dr 



(A5) 



where the matrix X is 



X 






1 





c 


B 


E 


A 




^A 











E 





C 


A 


^A 



(A6) 
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and the 2-vectors 



Kgfc I 



\ V M2 



kqQ 



(A7) 



with 



h" = 



M4 + g2g2 



- gQ) ( qQa + bM^ + 2ipM^a + _ ^2^2 



(A8) 



d = 



^o + (l + 2p)(^-l)^i, 



uj - qQ ' 



gQ(2 + n(n-3)(l + Q^)) 2a;a 4 + 

0^ 13 ^ uj-qQ ' 



a = 



n(n-l) + (n-3)(n + 2)Q2 



2/3 



^ = (n - 1) + (n - 3)Q2 . 
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